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do not rest on well-specif led sl:atlBtiaal models, Receiit 
approaches based on path analysis suffer from the defects of 
variable standardisation and lack of attention to measurement 
error , Ttim paper formulates a jnieasureinQnt model for quanti- 
tatli/ely scaled multivariate panel data. The model is applied 
to a data set indeKing two con struct s measured at three time 
points, mitlple sieasureraent of eaoh construct in corijunction 
^Ith tke meaaurament model allows eptlmatloii of a true v^ariance 
covariance maitrlx, Aiialysis of this matrlK producaa substantial 
ly different i^terpretatiQns of variable Influeiica than similar 
analyses of the original data^ 
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1 . Problems in the Inter pretation q£ Panel Data 



Sociological and social^psychological research 'workers have 
long been concerned with the attribution of causality to 
variables representing basic socio logical or psychological 
concepts. Given the difficulty off variable manipulation in 
many real social settings, sctentists have turnefl to statistic-= 
al methodology rather than to ejcperlmental techniques of 
investigation for help. 

From La^arsfeld's (1948) early dlscussloTi and exposition of 
the slKteen^fold table technique for qualitative data to 
Campbell *s (1963) and PelE and Andrews' (1964) simultaneous 
conception of cross-lagged panel correlation procedures for 
quantitative data ^ practicable data-^a^alytic technigues have 
been available for causal attributiDn* These technlgiieS/ how^ 
ever, have been hemvily criticised by statistical inethodo^ 
logists (e,g, Goldberger, 1971) anfl are consequently clouded 
with amblguitiea of interpretation!* One reason for the persist- 
ence of ambiguity is the lack of i^el l^def ined mathematical 
models to serve as foci for disGussiqn and bases for critique. 
Our Intention In this paper is to specify such modeLs as well 
as to eKposlte technique* 

In this paper v/e are focussing solely on qiiantltatlve variables. 
Consequently, the siKteen-f old table technique la not directly 
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relevant, Comparison of cross-lagged pan«?l corrcKitions w.Xf? 
the first widely advocated non-experimental technique for the 
attribution of cauaallti^ to quantitatively sc^aled variables 
In social-psychological research. In addition to the lack 
of a ciearly stateca statisticaL model, one universally re- 
cognised weakneas of this techriique li its serious distortton 
by comnonly occairring measurement erKOrs. One especially 
distorting event is systematic change in the reliability of 
variables over time. Such events occur, for example, wliea 
a true variable's variance changes, although the quality of 
the neasurenient remains the same. This is nearly always the 
case when any change In a true variable takes place (Wiley 
and Wiley, 1970). Therefore, it is difficult to justify the 
\ise of technique in typical social research settings. 

Path analysis is growingly accepted m a powerful frameTOrk 
for aisecting social data, it originated in genetic studies 
early in this century and has been systematically gaining 
favor among social scientists as a useful devise for strlpplag 
data to bare their (injcoiisls tencies with complex theoretical 
assuniptions and hypotheses. Path, analysis has the advantage 
of being completely specified in matheitiatical form. It Is, 
therefore, easy to criticize on both aubstantive and raetho- 
dologlcal grounds. Path analysis has the advantage of being 
multivariate and capable of variable aijgmentation while main- 
taining the basic asymmetries which are required in the 
simulation of causal as opposed to simply relational networks. 
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The appLicatioii of path analysis to quantitatively scaled 
panel data has been systematically developed by Duncan (1969), 

However, the path-'analytic framework^ as it is usually Imple-^ 
mented and specified in applied research work, suffers from 
at least two defects/ especially, when oriented towards longi^^ 
tTjdinal data^ standardization and measurement error. Standardia^ 
ation means rescaling of variables so that their standard 
deviations la the sample are in every case equal to one. This 
practice has bean generally criticized as wasting valuable 
information even in Gross-sectional data (Tukey^ 1954/ Blalockj 
1957a and b) , In the conteKt of longitudinaL data^ these 
difficulties hecoitie even greater, since most such analyses are 
oriented tov^ard the assessment of and the distinction between 
stability and change* Such standardization complicatea the 
detection of change or stability in the structure of the effects 
of causal fac±Drs# because invariance In the structure of the 
causal network^ which appears in the form of constant values 
for causal coef £ icients when the variables are not standardized * 
will result in non--constant values for those coefficients when 
they are standardized. 

This problem becomes even more difficult when errors contaminate 
the measurement of variables because there is no longer a 
simply specifiable relation between the two kinds of co*- 
efficients, rh^se ineasurement errors also have serious con-^ 
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Sequences for the estima tion of coef f iclDntti . Wi le^y and Wlloy 
(1970) and Wiley (1973) have discussed some of these distort-- 
ing influences within the conteKt of univariate loncfitudinal 
and multivariate data^ respectively* A general discussion of 
the complex consequences in a generalized regression frame- 
work has been given by Cochran (1968; 1972). ignoring measure- 
ment error in regression analysis complicates tba regression 
coefficient "attenuation" problem, commonly encQUntered in 
a simpler bi-^variate context, when measurement «arrors con-- 
taminate variable assesinient . Thio co^pleKif ication implies 
increases in the sizes of some coefficients as well as decreases 
(attenuations) in those of others. The compleKlty ot distortions 
in these key-indices^ used to make causal attributions^ effect-- 
ively distroys any hope of a simple rank-order ^relation be- 
tween observed and real coefficients. This implies that it 
is necessary to explicitly incorporate measurement errors in 
the formulation and specif ication of the model that serves as 
a basis fpr data-analytic procadures * When a single psychos- 
logical variable is errorfully measured this has been known 
as the problem of "the Jiiaasurement of change" in the psycho-' 
metric literature (see toonbach and Fur by ^ 1970) , This paper 
is one In a series which has the Intent of expoaiting models 
for the general analysis of quantitatively scalad data with 
measurement errors , 



A basic statement of tha fundamentals of such a model was 
given by Wiley and Wiley (1970) in a critique of a standardized 
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path analysis specification for unlvajriate longitudinal data 
(Heise, 1969). Wiley (1973) subsequently forniulated a class 
of such models for an econometr IdiLj^ly sophisticated audience, 
Keesling (1972) applied this baaia model to mill tlv^ariate cross^ 
sectional data and incorporated th^ notion of cgrrelations 
among errors. This was the first conjunction of Cronbach's 
multivariate generalizability thawy {Cronbach et al ^ , 1972) 
with complex multi^relational mod^l^/i common in econometrics. 

This paper presents procedures an<^ techniques based upon 
Btatistlcal models which incorporate more than one theoretical 
concept meaeured at more than on^ point in time . The model 
includes measurement errors in th# Otb^erved variables which 
may be correlated over time* The ^Kposition and diBCUssion 
is oriented toward empirical res#a;rQh workers in social and 
psychological science* We hope to :rtiake the models and procedures 
both compreh^nBlble and practicabX^ to research workers who 
have not received extensive advanced training in statistics, 

2 . A Basic Measurement Mod^l for O i jia.ntitatlve Variables 

In a data set representing two coj^atructSi each measured at 
two points in timer there are sev^aral relatively simple pieces 
of information which might be deiirM in attributing causal 
impacts. In the absence of measu^TOent error, there are three 
quantities which describe simple tftlitributional characteristics 
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Of the constructs at the first time point: the two variances 
and the covariance. There are four more characteristics 
which represent the productive relations between the constructs 
at the first time point and at the second i the two regression 
coefficients relating Construct 1 at time two to the two 
Construe La at time one and the two regression coefficients 
relating Construct 2 at time two to those at time one. Finally, 
there are three quantities remaining: the two variances and 
the covariance of the residual parts of Constructs 1 and 2 at 
time two which are not accounted for by their time one versions. 
These ten quantities are equivalent to, i.e. contain the same 
information as and only constitute a reorganization of, the 
ten origina^quantitias: four variances and six covariances. 
They eonstiWte the sunmary of the distributional character- 
istics of the four measurements! two variables at each of two 
time points. Thus, In the absence of measurement error, we 
seek ten new quantities from ten old quantities, a possible 
task. When we add raeasurenient error, we must add new quantities 
to account for the characteristics of those errors. Clearly, 
when we do this, we increase the number of desired quantities 
beyond ten and thus eKceed the number of available pieces of 
inf oritiation . 

Researchers have made a variety of assumptions so that the 
number of values they wish to estimate is no larger than the 
number they possess. The simplest such assumption is that 
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there is no irror of: measuremerit. ri all the etror %rarianceB 
are tmto ^ theii yrm onLy ^eed es-tlmate ani r^axrang^ the origin- 
al vaLues: the four" xrtriaraes aiid ttie slx co^aria^c^s, Therre 
is a aitect CDJrespon4eiice bet^^een the o^rlgimL valines calculat-- 
ed arid- the xrhliimm disired^ irnSortTOateLir , such an asi^ptlon 
is raraLy^ legitimate in BociaL icienceaiid anaLyses hased on 
it are offten mislia^diaig . 

A secoM approach has hm%n to Lnc^ease t3ie number of i terns of 
Information available beyond tl%a ^urciber of quantitiess one 
^litieB to assess. exaniple off tliis apparoach is Hei^e^e e%- 
paasioJi of t^^o^time point iiodeLs foir analysis oJ three— tlma 
point panel data (1S69)* H€ asames that all the E-alatioriahips 
artiong the vaxisblis may be e2<pL4i^ea postiila+ing thait a 
eimple Maarko^ nnoail haias, l^e, that v^ariatolea at a specific 
time are datermin€d only^ tho^e at the prscaddiig' poirmt. He 
also assuTOS stable reLtabilitics for the \rarlaMes, Wi^thcut 
going into datalL / this iKBansien gii/es Hilse enoixgh addition-^ 
al values (variances and cemrianca^! to Identify a roodei wWah 
incluaes ireasurement error - Whiae an improvement ovet one witht- 
out iiieasia:cimeiit error^ this inodel is etill ^ulte r estrictive. 
A dlseiaesion of the model' s atllity can toe found in Hii^e (196 9), 
and ctittcism of some of its assiiiriptioni in Wiley aiid- Wiley (1 970) , 

We bell we that tha^i is a more iffectim wsy tc^ increase the 
amount of infformtlo^n ieyorad theniimtoerof quantltiis one fteed^ 
to eetimti. Instead oi ejcpanding i:ha numler of tima joints^. 
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•we suggest expanding the number of measxires of each variable 
at each point. Two measures of each variable at each time 
point produce enough additional information to allow the 
calculation of all the quantities in a quite general two-time 
point moael. We can even calculate all of the important charact€ 
Istics of the true variables at a siiigle point in time. This 
will be Illustrated in the next section. In fact, ve have more 
Informatioii (observed variances and oovar lances ) than we need. 
This extra information also permits us to assess the adequaoy 
of the model. Increasing the number of measuies of each variable 
beyond two, permits further loosening of assumptions. For an 
extensive treatment of the multiple measurement of psycho- 
Logical constructs see Cronbach et al^ (1972). 

In order to specify the model, it is necessary to make some 
assumptions about the structure of the measured variables. 
Specifying these guides us in designing our measurement oper- 
ations and variable definitions so that they more closely con- 
form to these assurnptlons . 

A traditional way of making these specifications is, that each 
measurement be a linear function of the underlying (true) value 
and an error ^ and the errors be uncorrelated -with the true 
values for every measure and errors in other measures. This 
means that a covartance between two different measures of the 
same trait should resuit entirely from the common true variable. 
A covar lance between a measure of one trait and a measure of 
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another should only leiult from the true relation between 
the traits. ^- 

Iii describiiig the nteasurenient characteristics of oiir model, 
we will begin mth tw coiistructa measured Im t^o Yerslona 
at a single point in tiine. Denoting this time point by t and 
the version by "1% the model for a datum resulting Sxom the 
first operatloiiai definition of the first coristruct iii 

(2.1) Ki^ ^ 4t ^ ^^t ' 

where xi^ represents the obierved rtieaaurenient , 5^ the true 
value ^ and the error of measuremeiit , We assunte that tha 

error is uncorrelated with the true values of its own and 
other variables as well as vith other errors at the same 
time poiiit. Siniilarly, the model for such a datum, of the 
second CDnstriiat is : 

(2.2) yi^ ^ ^ ^^t 

where yi^ represents the observed measurement, the true 
value, and Si^ the error of measurement, EKamintng the vari- 
ational and co^^ar^ational structure of the obierved variables 
in this simple model, we obtain i 

(2.3) Var(X3. ) ^ o| + 

(2.4) Var(yi ) + cjj , and 

■t ^*t 

(2.5) Cov^(:ki^, yi.) ^ o 

^t,^t 
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for the variances and covari^nce of the observed variables. 
These three observable s axe cQmposed of fiv^e basic eleiments. 
The fact that three obserTables are composed of five elements 
implies the ImpossiblLity of the evaluatioii of those elements 
without additla^al ijif crmatton. 



rf we f however, take active trole in designing the data 
GolLecticn process |. we maj generate aflditio^nal Information 
^hiah mil alLov? us to evaluate these uiiobservabLe v^artanaes 
and cov^arlances. Wiley ( 1973> has suggested that soGiaL measure^ 
ments be designed according to strategies similar to those used 
in psychologtctl measurement ^ that is^ by genexating more than 
one operational version of each construct, in addition to (2.1) 
and (2.2) abovei^ we might assume that each of the construots^. 
at tiraa t, 5^ and rr^r has a second operational deftnitioni 
(2.6) - 6^1^ + e^^ , and ^ 

where and represent the possibilities of differencas In 
the scales of measuieiTient hmt-vmem the first ana second -versions 
of each coristruct at tlaie t. ^hase addi inal obsermble 
Yarlablas geneEate two more ofcsarvable Vi. icas and five 
more observable covariances : 



(2.6) Var(x2. ) * pfa| + o-^ 



(2.9) Vair(ya, ) =^ w^a^ + 
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(2.10) CovUi^, 5ci^) = e^o| 



(2.11) Cov(3{a, , yi.) « S.Cp ^ 

(2.12) Cov(y2, , XI ) = a,o- 



(2.. 13) Coviyz^, Yi^) = ot^o^ , and 

t 



(2.14) Cov(k2.^ Y2^) = a.B^fft ^ 
t t t t i^/Ti^ 



lasgectlng the rlght-hana aids of theee ten coinpOBitlcnal 
equations, wm find niae unobseE\rabla coinpDiientsi ai , ^ i 



c^L f ^ ^ S ; ^ and o| , As Long the 

^ife ^t ^ ^^t -^t 



guantitieE ana whLch are coef Eiclents characterizing 
th^ raeasuramerit scales of the ^?arlablai/ aie iiot zero (a Eexo 
value implies that the new measures do not realLy relate to 
the appropriate constructs) aiia as long as the two coristriiGts 
are related (cr is not equal to sero) , this system of ten 

equatdons with nine unknoi^ns may be solved and the unobservable 
Qompoxiente becorae aalcula table from the observable ones. 

This ahaiige in nieasureinant design has allowed us to separate 
the interrelations among the observed constructs into those 
three aapects specifically d^e to the underlying CQnetructs 
ana those six due the dlstortioais of meaeurement errors . 

o ' m 17 
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2. 1 Example 

In this section we will present an example vhlch will illustrate 
a procedure for estimating the variances aind .covariajiaes of 
true and error variables. The data come from a study of 
conmunlcatlon processes oonduotea in Central America, fhm two 
constructs, we have chosen for Illustration, are televiilon 
watching by children and television possession by their fawillei. 
Each construct was measured twice at each of three points in 
time. The saine four questionnaire items were used at each time 
point. Coding of the responses to each Item at each poljit 
in time resulted in twel-ve variables. The covariance matrix 
for these variables is displayed in Table 1 . 

Insert rabli 1 abo-ut here 

We will demonstrate the estimation procedure based on the moael 
for a sing-le time point, presented above, using the data for 
Time Point i. Since the number of quantities to be estimated, 
nine, is fewer than the nunber of observed guantitlei , ten, 
there are several possible procedures. Two ooinnioiily chosen 
methods of "optimal" estimation are generalized least squares 

mascimun likelihood (see e.g. Goldberger aiid Joreskog, 1972). 
We have chosen a simpler method, for illustrative purposes, 
which consists of eliminating one of the eguations ana solving 
those remaining. This method produces a relatively but aiot 
fully affiaient estlinate acoording- to the above criteria. 

18 
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The renaiiiiag equation can be used for the purpose of testing 
the adeguacy of the model and we x^lll illastrate this use 
aLso* 

Since Time Point 1 w^as ohoBmn for tllustratlon ^ tlie first four 
hy four submatri^c in Table I is the basis for cur calculations-. 
2rk ail that foLlcws, k wlLl characterize obm&tvmd television 
pDSEeision while ^ will denote observed telev^iatoii wmtcliljig. 
Recalllag from above that the first siabsaript specifies the 
verston of the construct vrhile the second represents the tiine 
point at wh-ich it is meaeared^ the variables^ In crder^ are 
symbollzeai ^ti r Xa i , Yitr and yti, which represent posses alon 
v^ersion one (time one) ^ pDsseiSion version two (time anm) , 
watch ing verston one (time one) / and TOtohlng x/ersloin tv;o 
Ctline one) , respectively^ 

According to Equation (2.5)^ 5^ ^ ^ covCxat, which froiri 

rabie 1 Is egual to ai9.''^ Jron Equation C2,L) , $i ^ Bih / 
^f^n, ^ oav(Hzi, yii)/Q' which is equal to (0 ,443/0. 479^^) 

0,924 v^hich is the metric-coefficient for the second indeK 
of televasicn possession in the scale of measurement of the 
first. Since both of the indices were hased on dlchotoinoiis 
Items # where the natural interpretation of the positive 
alternative for both items v?as possession , the scales of 
measiireinent should be almost the same. If thie were In fact 
the case^ 'pi would be egual to 1*000, The value which tve Jiave 
obtained is therefore quite reasonable, Slimilarly, Srom Equation 
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(2.12), Si = ai^|,_n/^Ciji3i ' °Sv(y2i, xn)/6^.^^ which is 
equaL to (O. 312/0. 479«) 0.651, the analogous coefficient 
for television watching. The fact that the estimate is less 
than one is also reaspiiable since the range of varimimi of 
the second vatchlriq item is less that of the first, yniitinu- 
ing in like fashion, fyom Equation (2.10), 6^ = SiO^ /Bi 
= cov(5<2i^ xit)/Bi equals (0,227/0.924=) 0.246, the true score 
variance of talevislon possission. For television watching, 
from Eguation (2.13), ^ « a^o'^^/a,^ c6v(y2i, l/Si equals 
(1,225/0.651=) 1.882. 



For teLe-vlslon possession the variances of the errors of 
itieasurement may be Qstimated using Eguations (2.3) and (2.8) 




(0.247 - 0.246=) O.QOlr and &^ = (f}o| 

£21 f I 



= var(X2i) - (ei)^8|^ equals C. 248 minus (0 ,353) (0, 246) which, 

f 

in turn, equals 0.038, Theie small values indicate that tele- 
vision possession is qiiita reliably reported by the children. 
The variances off the errors of measurement for television 
watching may be estimated in a similar fashion using Eguations 
(2.4) and (2,9) as follows: 




equals (2.313 - l.€S"2«) 0.431, and 

^'Ui " ^^^T"^^^ - ' vSr(yai) - (S,)^&^^ 

equals 2.264 minus (O. 423 ) (1 882 ) which, in turn, equals 1.468. 
These valiias imply that the measures of television watching 
are considerably more eerorful than those of television possession 
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even considertng the differences in measurement metric. It is 
also clear that the second index of watching , which is based 
upon long term frequencies of watchtiig particular television 
prograins, is much more error ful than the first, which asks 
about general watching in the pracedimg veek. 

Recalling that we have not used Equation (2*14)^ covCKai; ^gj) 
= aipig^^^^^^ for estimation purposes, we may ennploy it to 
assess the adequacy of the model, ThLm may be ace oinpli shed 
by comparing the tabulated value ^Ith a prediction baaed on 
our estimates of its components, The value for this CQvarianoe 
from Table 1 is 0,288, Predicting it Crcm the estlrnateB ai, , 
^iifni' obtain (O, 65l) (0, 924 ) Co ,47 5) or 0.286. Slrice the two 
quantities are alinost egualv we a^e pleased with our assess-- 
ment of the adequacy of the modal. Table 2 reports, for Time 
Point 1, all of the estimates and the ratios of the observed 
versus predicted covariances. 

Insert Table 2 about here 

We may, in the same way, estirnate the corresponding- quantities 
character Izlng each of the other times of measurement. These 
are based on the central and final tanx by four subcovariance 
matrices of Table 1 , The results Sox Time Point 2 and 3 are 
also suimnarized In Table 2, In each case, the covariance ratio 
ijidicates that the rnodel is quite adaguate for our data. 
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It shouad be noted that the model, as specified, assumes con- 
stancy over time in the metric off the first version of each 
variable. Since we do not include parameters such as a and g 
in these versions, we make no allowance for scale changes in 
them. We do, however^ allow such changes In metric for the 
second version of each variable. Accordingly, we have sub- 
scripted the a and p values separately for each time point. 
Ibis provides for a different discrepancy between version one 
and version two to obtalri at each timm point. Consequently, 
general metric inconsistencies ffor thia second version of each, 
variable are possible. 



This is logically necessary for the watching variable. The 
first version of the watching variable has precisely the same 
operational definition for each time point. However, the tele- 
vtslon prograiTiB, used to specify the Second version, change 
from time to time. Becaiise of these changes, the metric co- 
efficients should vary. On the other hand, the possession 
variabLe has constancy of definition £or both of its versions 
at each point in time. As a conseguence, we would not expect 
tlie true values of these metric coefficients to vary over tine. 
If we inspect the estimatas in Table 2:^ we, in fact, find 
slightly less variation in the coefficients for possession . 

We can also assesa the precisions of the nieasureinents and the 
wajfs In which these precisions change over time. The reliability 
coafficlent of a measure Is usxially dtfined as the proportion 
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of variarice attributable to the true variable. Table 2 conveyed 
the variances of the errors for each meas\irement operation at 
each time point. We may use these to aasesa the reliability 
of each rneasuie since we also have estimates of the variances 
of the true variables and of the scale parameters, 

Wiley and Wiley (1970) have strongly argued that this ordinary 
coefficient Is mlsleaaing. The standard error of measurement 
(square root^ of the error variance) defines more adequately 
the a^uracy of a measurement since it inaicates, in an average 
sense, how far an ohserved measurement is likely to be from 
the true one. Since the standard error of measurement is not 
the only characteristic which influences the reliability^ the 
traditional coefficlant can vary because of changes in the true 
status of the yariable as well as changes in the characterlstios 
of the measurement operation* Concretely^ as the variance of 
the true variable increases, other things being equals the 
reliability increases. 

Table 3 displays the reliabilities for each measurement at each 
point in time. The measurement properties of the variables, as 
indexed in this traditional way, vary widely. Since the true 
variances of the possession and watching constructs change over 
time, especially the latter, we find f luctiaations in reliabilities 
which are not only due to changes in the average si^es of the 
errors but also to changes in the distributions of true watchlncr 
behavior on the part of the children. 
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Insert Tatole 3 about here 



3. A Specification of th$ Measuremen t Model, ftllowlng 
Cgrrelated Errors for Repeated Measurem «3nt3 

Just as we developed representations of each of the observed 
variances and covariancea for each time point, w© can construct 
similar representations for the cross-covarianoea between time 
points. Before, we were interested in defining ^ tm example, 
the covariance of j and ya j , both jneasured at fim© Point 2 . 
Now, we are interested In defining, for example, the covariance 
of Xj i and Yzz> the first measured at Time Point 1 and the 
second at Time Point 2. 

Vfe specified that the errors were uncorrelated with true 
values and with errors in other variables. We did not specify 
the relations among the ©errors in the same variaWe at different 
points in time. Since it ts reasonable to expect that the error 
in the measurement operation for a specific conatruot at one 
point m time may be correlated with the error Iw that same 
operation at a later point, it would seem advantageous to 
allow for this in the spe«lf ication of the modeX. More concretely. 
It seems reasonable to aa»ume that a child who owrgtates his 
general television watchiag at Time Point 1, may «l»o overstate 
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that at time two* It is possible to allow this aaditional 
flexibility because there are many more observed variances 
and aovariances than values which 'v^^ desire to estimate. 



Of course^ the allowance of correi^tiad errors over time for 
specific measurement operations impll&s that the covariances 
among the resulting measurements ^111 not directly reflect 
true score covariances as they do within time periods. For 
example: 

(3*1) C0V{Kll rXl2) 

^ Covdi + eii.la + eia) 

^ Cov(Ci,5a) + Cov(|i,eiz) ^ Co^Ceai.Ca) + Cov(ei|,£ia) 
CifSa en, Eli 

Vfith this new specif Ication in mindi/, we can turn to the general 
task of detailing the forty^elght ^yrtrtsolic compositions of the 
cross^covariances between time points (alKteen Cor each time^ 
point pair) * These are summarized together with the thirty 
previously specified within^time poli^t variances and covariances 
(ten for each time^point) in a symbolic varlance^covarlance 
matrix (Table 4) * 

Insert Table 4 about h©r^ 

We have forty-^lght new pieces of iinifiormatlon ^ but desire only 
twenty-four additional quantities^ thi-f twelve crose^covar lances 
of the true vaarlables (4 for each of 3 pairs of time points) 
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and twelve error cross--covariances (also 4 for each of 3 p.riirs 
of time points) * 

In Table 4 the symbolic representations of the crosi-GovarianceB 
of the observable variables are located in the three off-^diagon^ 
al blocki . The diagonal blocks symbolically represent the 
variances and covarianGei within the three time points. This 
sexpartite structure parallels that of the empirical values 
for the varianaee and covariances displayed in Table 1 * With 
the thrse-^time point data^ the number of available pieces of 
Infornnation exceeds, by more, the number of desired quantities, 
than in slngle-^tlme point data. This Implies that there are 
many more ways of estimating these quantities. A statistically 
more elaborated solution to the problem of estimation than 
that given below would make more efficient use of this large 
amount of additional information 

3 , I Example 

This subsection illustrates the computation of estimates of 
the aross-covariances for both the true variables and errors* 
We demonstrate these computations using the cross-covariances 
between Time Point 1 and Time Point 2. We will refer to the 
entries In Table 4 by the row (1 ^ 12) and column (I ^ XII) 
numbers corresponding to the particular cross-Govariance . 
Since coitibinatlons of Time Points 1 and 2 are in the first 
off -diagonal block, the row numbers (Time Point 2) will always 
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be between 5 and 8, while the column nuniberB (Time Point 1) 
will always be between I and IV, 



Two off ^diagonal cells (6, 1 and 5^ II) in the off-diagonal 

block include the cross-^time covariance of true possession , 

0^ ^ , Referring to these two cells in both Table 4 and 

Table 1^ we find that BiO^ ^ correspondi to 0.201 and B%a^ ^ 

corresponds to 0.196. Since we can obtain the estimates of 

Ba and 0i from Table 2 (0.998 and 0.924), we may estimate ^ 

by 0.201/0.998 = 0.201, or by 0.196/0.924 ^ 0.212. For the 

final estimate, we choose the mean of these values: 6^ 

equals 0.207. 

The diagonal elements of the off-diagonal block contain com^ 
ponents representing both true and error cross^time covariances 
For example, &^ , ^ ^ ^ + a = 0.198 (5, I). Sub- 

tracting the covariance of the true variables, estimated above, 
we obtain 9 = 0.198 0.207 = -0.009. Similarly, 9 

^ Bi&iO^ p ^ = 0.205 (6, II), therefore, S 

^ 0.205 - (0,924) (0.998) (0.207) ^ 0.014. Applying the procedure 
for estimating ^ to the entries in cells 7, IV and 8, III, 

we obtain for watching g ^ 1.602. Consequently, 9. . 

- -O.I43 and 9^ ^ ^ 0.617. 

0 2 1,022 

If we inspect the eight remaining cells in the. of f =-diagonal 
block, we find that half of them can be used, in an analogous 
fashion, to estimate and the other half to estimate 
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The latter is a component of the following cells: 
Cell 7 , I I ^ = 0.402 

Cell 7, III eTa. ^ ^ 0.375 

Cell 8, I I a^Op ^ = 0.212 

Cell 8, III fioia^ ^ ^ 0.185. 

Dividing by the previous estimates (0i = 0,924 and Sa 0.571)^ 

we obtain four distinct values for i 

0.402 0.402) 

0.406 (= 0.375/0.924) , 

0*371 0.212/0.571), and 

0.351 {= 0.185/(0.924) (0,571)) . 

The average of these four values 0.383 ^ 8^ . Followinq 
the same procedure with the remaining cells we obtain 0.441 

62 f ni 



We have now estimated the four true crosa-»Govariances relating 
possession and watching between Time Points 1 and 2. We have 
also estimated the repeated--measures error^-oovariances for the 
same two time points . The procedures we have used to produce 
these estimates can be equally well applied to the cross-co- 
variances of the other two pairs of time points. The resulting 
error covariances are displayed in Table 5. If we combine the 



Insert Table 5 about here 



estimates of the true cross-covariances for each pair of time 
points with the estimates of the true variances and covariances 
within each time point, we may compose the total variance-co- 
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variance matrix of the six true measures of the two constructs 
at the three points in time. This estimated six by six co^ 
variance matrix is displayed in Table 6* 

Insert Table 6 about here 

This matrix serves as a new basic data set for all subsequent 
analyses of the true measures. When we compare the results of 
analyses based on observed data with those based on true 
measures/ we will first analyse the observed values in Table 1 
and then perform a comparative analysis of the eBtlmated true 
values in Table 6. 

4 • Multi-relational Models for Longitudinal Data: 
The Effects of Measurement Error 

The variances and covariances of Table 6 allow the computation 

of unstandardiEed regression analyses relating the true versions 

of the possession and watching variables over time. We, in fact, 

have all the information needed to perform multiple regressions 

relating any variable to any selection of other variables. The 

only difference between our regressions and ordinary ones which 

are based upon individually observed values, is that we perform 

our computations from summary characteristics (variances and 
9) 

covariances) and that those summary characteristics are 
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estimated only indirectly from the original data rather than 
computed directly from the observations themselves. 

Regression analyses, analogous to those interrelating the 
true variables, will also be performed using the fall-ih] ^. version 
of those variables directly resulting from the original measure- 
ments* In computing the regression analyses involving these 
observed variables, we will use the variances and covariances 
of Table 1, These values are also summary. They were computed 
directly from the original individually observed values. 

Another distinction between all of our regressions and those 
usually performed with individual values is our omiBSion of 
a constant or intercept term from the specification of the 
regression model. This omission has no important conBequences 
because the computation of regression coefficients from variances 
and covariances automatically eliminates the effects of variatloni 
in the means of the variables. The regression coefficients 
remain the same as they would have if the constant term had 
been included. The analyses will relate both television 
possession and television watching at adjacent time points. 

We will relate television possession at Time Point 2 to possession 
and watching of television at the first time point. Substantively, 
we would eKpect a close relation to SKist between the two 
measures of television possession, i,e* we do not eKpect families 
to dispose of their television sets. We might expect a small 
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relation to ejcist betwfsen prior watching and later possession 
of television since those families who do not already own sets 
but have children who watch television outside of the home^ may 
over time be influenced to acquire a set. 

The additive regression model, in this case, may not be most 
appropriate^ because we do not ejcpect a relation between watching 
and subsequent possession for families who already possess tele= 
vision sets. I.e. our illustration of the methodology is only 
approximate with these data since the possession variable is 
a true dichotomy. 

If we examine our eKpectations for watching , we would expect, 
with or without prior possession , a close relation between 
television watching for adjacent time points. For two children 
who watch television equally long at the first time point, we 
see several possibilities for the relation between possession 
and subsequent watching . It is, perhaps, most likely that any 
effect of television acquisition at Time Point 1 would have 
Immediate consequences for television watching. Therefore, all 
effects on subsequent watching would be mediated through 
watching at Time Point 1 and, once we have taken this into 
account, we would expect no relation between television possession 
at that time point and later watching. However, if we had some 
speculation that the full effect of television acquisition on 
watching behavior were not felt iimnediately , but only later, 
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we might expect an independent influence of possession on 
subsequent television watching even when we allow for prior 
watching , 

We may test our expectations by first performing regression 
analyses relating the Time Point 2 values to the Time Point 1 
values using the observad measurements* That is, wa proceed 
in a fashion that is usual when such regressions are carried 
out without adjustment for unreliability. This implies that 
we fit the models i 

(4a) Ki2 ^ ^13,^11 + '^i^yii + Eia^ , and 

(4.2) yia a Yl,.xn + llVil + 

X y y 

in place of the models; 

(4.3) fa = Xgli + A^ni + 8 , and 

(4.4) na = Y|Cl + Yj^ni + 0 

and treat our estimates of Xi__, Xi^, yj_^, and yi^ as if they 
were estimates of X^, y., and ■ Using the values in 

Table 1 , we obtain the following estimateB ; 

possession; = 0.730 (Possession), ti^ = 0.046 (Watching); 

Watching; = 0.665 (Possession), * 0.494 (Watching). 



If we compare these estimates with our prior expectations under 
the assumption that these numbers represent the relations among 
the true variables, then we see that our expeotations for 
possession are verified. There is a relatively large impact 
on television possession at Time Point 2 of that at Time Point 1 
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and there is almost no relation betv^een watching and subsequent 
possession of a television set. 

The results^ hav^ever, for television watching at Time Point 2 
were less expected* Although we do find a moderate effect of 
earlijr watching , we discover a close relation between initial 
possession and subsequent watching of television. This seeming 
"sleeper effect" of delayed influence was only speculated about 
before, but seems quite apparent in the data. 



Let us now compare the results of these fallible regressions with 
some based upon the estimated interrelations of the true variables, 
The estimated coefficients are: 

Possesslon i = 0-735 (Possession), %^ = 0.042 (Watching); 
Watching i ^-0,004 (Possession), = 0,841 (Watching), 

Inspecting the possessions regression, we find almost no 
difference between our new estimates and the old ones, which 
indicated that watching had little or no effect on subsequent 
possession of a television set* However, whan we look at the 
determinants of watching , we find a striking disparity* There 
is no relation at all of possession to subsequent television 
watching apparent in our new estimates, and the size of the 
coefficient relating initial and later television watching 
has almost doubled. The natural interpretation of these equations 
is, in fact, that there are no causal inter connections between 
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television watching and possassion. There is no surprising 
"sleeper effect". 

If we eKamine the discrepancies between observed and true 
relations carefully^ we may divine the cause. The accuracy 
of measurement of the possession variable at Time Point 1 is 
quite high. This may be seen from the reliabilities in Table 3 
or the variances of the measurement errors in Table 2. This is 
not true, however^ of the watching variable. The variances of 
the errors of measurement are substantively larger than those 
of the possession measurements. Even taking differences in 
metric Into account^ the reliabilities are substantially lower. 
The effects of these errors # on our assessment of the determinant 
of possession at Time Point 2, are small. The large effect is 
that of the most reliable variable', initial possession ■ The 
small effect is that of the unreliable variable t initial watching 
As a consequence, the larger errors in the initial watching 
measure:T:ent have little impact on the regression because the 
true watching variable, in fact, influences possession very 
little. 

The case is considerably different, however, when we look at 
the regression coefficients for the determination of watching . 
Here, the most precisely measured determinant has no influence, 
while the least acourate eKplanatory variable has a very large 
impact. Since the true explanatory variables are positively 
related, error of measurement in the most powerful determinant 



(watching) not only attenuates the estimate of its effect but 
also weakens the "control" exerted On the relation between 
initial poss^^sslon and subsequent watching of television. The 
consequence is what might be called '-undercontrol " , A true 
effect of prior watching on subsequent television watching 
is spuriously attributed to possession of a set, because of 
the high reliability of the prior possession measure, the high 
relation between true prior possession and true prior watching, 
and the lower reliability of the prior watching measure. 

It should also be noted that in this case the interrelations 
among the errors in particular variables at the two points in 
time did not have a dramatic influence. These relations were 
small enough so that they had no material effect. If, on the 
other hand, they had been strong and positive, they would have 
partially compensated for some of the distortions In our 
estimated effects. However, this is only true, when there is 
no Influence across variables. In the case of true cross-varlabl 
iiifluence, large positive correlations among the errors in 
repeated measurements would result in too little attribution 
of influence to the other variable. 

We may replicate these analyses for Time Point 2 and Time 
Point 3 measurements. Again, we initially examine the observed 
regression using the first version of each construct i 

Possessloni - 0,792 (Possession), ti. - 0,011 (Watching); 

Watching I Yi^^ - 0.794 (Possession), - 0.467 (watching). 
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The pattern of results is obviously similar to that of the 
earlier observed regression. Turning to the true regression^ 
we obtaini 

Possession s = 0*917 (Possession), ^ -0,025 (Watching); 
Watching ; = 0.087 (Possession), = 0.857 (Watching). 

The pattern of change and the resulting implications for the 
interpretation of influence are the sanie as those we found 
above. It is interesting to note that not only do the true 
effects have the same pattern for each pair of adjacent time 
points / but also/ that the biasing effects of the measurement 
errors are the same. The replication did not help us see 
through the fog of error. Only an explicit assessment and 
accounting for those errors in our model and in our analyses 
helped us blow away the fog, 

5 , Summary CommentB 

We have reviewed some earlier attempts to deal with the 
analysis of panel data which have treated the problem of 
causality attribution to quantitatively scaled variables. 
Until recently^ most such discussions of panel data have 
focussed on procedures comparing cross^lagged correlations. 
These practices suffer from many defects/ whose severities 
have been difficult to assess in actual cam&B, because the 
procedures do not rest on well-^specif led mathematical or 
statistical models. 36 



Recent concern for more mathematically explicit statement of 
the bases for such procedures, arising out of genetj.cs and 
economics, has led to the increasing use of multi-relational 
models for the analysis of panel data. Also, the longstanding 
emphasis in individual psychology on error of measurement 
has gradually becoma more sophisticated* A combination of these 
concerns seems to be in order* The conjunction of measurement 
models, relating true variables to their measured manifestations, 
with multi-'relational models, relating multiple measurements 
of many true variables, is greatly needed. Measurement errors, 
have large distorting effects on the assessment of variable 
influences when they are not explicitly taken into account 
within a multi-"relational setting. 

This paper has formulated an eKplicit measurement model for use 
with the analysis of guantitatively scaled panel data. The 
model incorporates more than one measurement of each construct 
under investigation. This multiple measurement allows the 
assessment of the interrelations among true variables at a 
particular point in time. The addition of interrelations among 
the errors of measurement over time allows assessment of the 
true cross^time relations of the underlying constructs as well 
as assessment of the accuracy of the measurement procedures. 
Once the intra^ and inter--time point relations among the 
constructs have been assessed, multi--relatlonal models in- 
corporating them may be implemented. The resulting estimates 

37 



serve to assess the relative degrees of inter-construct inipact 
through time. This precipitation of true from observable, 
this purification from error, allows us to eliminate the dis- 
tortions which these errors cause when they are not taken 
into account. 

We illustrated these models and their implementation as practic- 
al data-analytic procedures by analyzing a three-time point data 
set with double measurement of each of two constructs. Our 
initial example allowed us to estimate the variances and co- 
variances of the constructs separately at each time point with 
no contaminating effects of error. Our additional specifications 
of time dependencies in the errors allowed us to estimate the 
cross-time relations among those constructs, as well. The full 
set of derived interrelations among the constructs at the various 
time points allowed us to fit multiple regression models to 
assess the true impacts of the underlying constructs on one 
another between adjacent time points. We demonstrated that 
there were large effects which had severely distorting influences 
on substantive Interpretations, These distortions were effectively 
removed^ however, by the procedures. 

Although we illustrated the models and procedures with data 
involving only two constructs and three time points, the 
general strategy is valid for any number of constructs and 
any number of time points. It should serve as a valuable tool 
for eliminating the distorting effects of measurement error in the 
analysis of quantitatively scaled panel data of much greater 
complexity. 38 
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Footnotes 



1) 



2) 



3) 



4) 



5) 



A dlBCUBsion of the issues basic to our thesis in the 
context of qualitative data may be found in Murray (1971) 



Campbell (1963, p. 240) has stated this problem clearly i "A 
variable which Increases in reliability from Time 1 to Time 
2 will, ceteris paribus , show up as an ' effect • rather than 
a ' cause' * " 



As stated in Wiley and Wiley (1970, p. 116) i ''Each atan^ 
dardized parameter is a function of more than one unstan- 
dardized parameter. In general, if two or more of the un- 
standardiied parameters of a model are eguai, the correspond- 
ing standarai^ed parameters will be unequal because they 
are not related to the unstandardized parameters by a n 
equivalent transformation." ^ 



While covariances between distinct measures are not affected 
by errors of measurement, under these assumptions, correl-- 
ations are attenuated, since they are defined by dividing 
the covariance by the product of the standard deviations. 
These standard deviations are inflated by the measurement 
error, and therefore, the correlations are deflated. 



Following are the questionnaire items, translated from 
Spanish with the coding of responses for each variable. 

Television Watching 

A, How many times did you see television in the last week? 

Coding 

none O 

one or two times 1 

^ three or four times 2 

five or six times 3 

every day 4 
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B, With what frequency did you see each of the following 
programs? 



Every Once or twice 

week per month Rarely Ne\ 

A, Tom Jones 

B* Tarian ~^ ~ ^ — 

C. The Office _ ~ ^ — ~ 

D* Bonanza " ^ — — ^ ^ 

C. Land of ' 

Giants 



Coding 



Codes for each program were suimned and the total was re- 
coded as follows s 

Sum Recode 

O O 

1.2,3 1 

4.5,6 2 

7.8,9 3 

10,11,12 4 

13,14 5 

16 6 

At each time point, program names were changed to corres 
pond to common preferences. This resulted in substantial 
differences in means and variances among the time points 

Television Possession 

A. Of the following Information media, which do you have 
in your home? 

newspapers 

magazines 

radio 
television 

^ L_ books ^ 



Children who answered "television" received a score of "J 
others received "O", 

B. In what places do you see television programs? 

in your house 

in a friend's house 

in a relative •s house 

elsewhere 4^ 
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Children who answered "In your house" received a score of 
"1", others received "0". 



Within the context of data from a single point in time/ the 
model may be "optimally" estimated by either generalized least 
squares or maximum likelihood, within the frameworks of Wiley, 
Schmidt, and Bramble (1973) or Jttreskog (1970) , A computer 
program for carrying out the analysis is also available 
(jareskog, van Thillo, and Gruvaeus, 1971) , The main conBequenc 
of such "optimal" estimates is that they are generally more 
precise (have smaller standard errors) than the more easily 
computed ones proposed here* As the models become more complex^ 
current "state-of --the^art" methods become inadequate* The ex-^ 
tensions discussed below, for example, cannot be "optimally" 
estimated using existing computer programs* 



The symbol """" refers to the estimate of the quantity over 

which it appears* When one such symbol appears over more 

than one quantity, the result refers to a single estimate 

of the composite quantity* When more than one of these 

symbols appear in a single expression, the result refers 

to the composite of the estimates of the individual quantities* 



See Footnote 6) * 



All of our regression computations were performed either by 
hand or using standard regression analysis computer programs. 
Many such programs accept summary as well as individual data 
as input. The summary input may take the form of variances 
and covariances or of correlations and standard deviations ^ 
each of which is easily convertible to the other . Programs 
which require means as well may be given arbitrary values 
in place of them and the resulting constant terms may be 
ignored . 
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Tabli I Cmtim Mmix for Two Vanions of IMm ?mmm and Wiiciig it Eich of Hiree Time Poiflii 
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Table 2. Estimates of Parameters for the TTiree Time Points 



Tuni Pdnt 

1 _2_ 



g .924 .998 .978 

S .651 .571 ,635 

^ .246 ,232 .238 

1.882 1.800 2.088 

0 

0. 479 .476 .502 

^ s n 

0^ .001 .018 ,012 

El 

.038 .019 .022 

of 431 .653 .387 



of L468 1.163 1439 

observed cov(x 2 * y 2) cov(x 2 ? y 2) 



predicted cov(x y 2) 



1.01 LOO LOO 



Table 3. The ReUabUiitiei of the Variables for the Three Time Points 

Time Point 

Construct Vernon Symbol 1 2 
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.928 


.952 
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.912 
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1 
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.335 
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Tabk 4. SyraboEc Covarianca Sttuctufg for Two Versioni of iM tmmHi 
at Each of TTuiaTima Points 
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T8bIe4contiiiued 



(1) 



50 



{3} 
(4) 

(6) 

(3) 
(9) 
(10) 

(11) 
(12) 



h ' C2 



0! + 0^- 

k hi 



0f f + 0. . 
92 43 El 2 'Ell 



I1 ' ^1 



50. 



ilpil2 



S2 



2' '12 



l2'n2 



^2' '13 



52*^3 



i 
I2 



12' tii 



a + 0, I 



Hi 022 



0,0 



a,0 



(2) 
(3) 
(4) 
P) 

(7) 
(8) 

(10) 
(11) 
(12) 



k hi 

hi 



^3 ni 



to 



1) 



%\t i Estiniated Vaiianc^cQvarianca Matrices of the MeESuratneni Eirors 
(CDfrelationi in Paientes) 



ConsMt, VoFSon Sy mbQl Time Point 1 
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TMiilue for the eofrelatlon k geater than one bicauie mpling vifiition 
has produced an abnonnally stniU value for the Time Point 1 wror ^liaiics. 
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Table 6. Estimated Variance-covaiiajice Matrix of True Televdsioii 
Possession and ^aichiiig at Each of Three Time Points 
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